Onset of a Propagating Self-Sustained Spin Reversal Front in a Magnetic System 
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The energy released in a magnetic material by reversing spins as they relax toward equilibrium 
can lead to a dynamical instability that ignites self-sustained rapid relaxation along a deflagration 
front that propagates at a constant subsonic speed. Using a trigger heat pulse and transverse 
and longitudinal magnetic fields, we investigate and control the crossover between thermally driven 
magnetic relaxation and magnetic deflagration in single crystals of Mni2-acetate. 

PACS numbers: 75.50.Xx, 82.33.Vx, 75.40.Gb, 75.60.Jk 



Raising the temperature of a flammable substance can 
ignite combustion, an exothermic reaction between a sub- 
stance and an oxidizer that results in a chemically mod- 
ified substance pQ. Deflagration is a self-sustained com- 
bustion that propagates at subsonic speed via thermal 
conduction; locally burning substance increases the tem- 
perature of adjacent unburnt substance and ignites it. 
Deflagration is a dynamic instability of a system gov- 
erned by local reactions and diffusion. These systems 
are ubiquitous in nature-from cell growth to epidemics. 
The study of these nonlinear dynamical systems reveals 
rich phenomenology-including traveling waves, dissipa- 
tive solitons [2] , and self-organized patterns [3j H] . 

Magnetic relaxation is a diffusion-reaction process that 
can develop an instability and proceed as a magnetic de- 
flagration. It can occur in a magnetic material prepared 
in a metastable spin configuration; here the reaction is 
the reversal of spins that release Zeeman energy and the 
diffusion serves to transmit the energy to adjacent mate- 
rial. The phenomenon has been observed as fast magneti- 
zation jumps in molecular magnets [5j[6], manganites [7] 
and intermetallic compounds [8]. In these systems, spins 
that are in a metastable state at low temperature slowly 
relax to a lower energy state releasing heat. In some 
circumstances, the heat cannot be compensated by ther- 
mal diffusion, resulting in an instability that gives rise to 
a front of rapidly reversing spins traveling through the 
sample at constant speed. 

In this Letter we study the instability that leads to 
magnetic deflagration in a thermally driven Mni2-Ac 
crystal. We explore how to control the crossover between 
slow magnetic relaxation and rapid, self-sustained mag- 
netic deflagration using magnetic fields. Based on earlier 
work [9j [10] , we propose a simple model and perform nu- 
merical simulations that provide a good description of 
the experimental findings. 

Mni 2 -Ac molecules are magnetically bistable [11] due 
to their high S = 10 spin and a magnetic anisotropy that 
provides a large energy barrier between spin-up and spin- 
down states (see Fig.[]Ji). In the presence of an applied 
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FIG. 1: a) The potential energy with discrete energy levels of 
Mni2-Ac spin Hamiltonian, Eq.JT] b) Schematic of the sample, 
the Hall sensors, the heater, and the directions of the applied 
magnetic fields. 



field H = (H±,H Z ) the simplest spin Hamiltonian has 
the form: 



-DS Z - 9I1bH z S z - gn B H^S^. 



(1) 



The first term defines the anisotropy barrier (DS 2 ~ 65 
K) and the second and third terms are the Zeeman en- 
ergy corresponding to the field applied parallel and per- 
pendicular to the anisotropy axis. As shown in Fig.[]Ji, 
the energy barrier between the metastable and the sta- 
ble state is the activation energy, C7, and the energy 
released by the reversing spins is the Zeeman energy, 
AE = 2gfiBH z S z . The transverse field H± mixes eigen- 
states of S z reducing the activation energy and leading 
to relaxation by quantum tunneling of magnetization, 
particularly at specific resonant values of magnetic field 
(H z = kD/(gjiB) — 0.45/c T, for integer /c), correspond- 
ing to the crossings of energy levels with opposite spin- 
projections. The longitudinal field H z biases the poten- 
tial well also reducing the activation energy and increas- 
ing the relaxation rate. 

The abrupt complete reversal of the magnetic moment 
of Mni 2 -Ac crystals was first reported by Paulsen and 
Park [12 and later shown to take the form of a propa- 
gating spin reversal front [5 . Magnetic deflagration has 
been studied by sweeping the magnetic field or by raising 
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the crystal's temperature. Thresholds for both tempera- 
ture [13 and applied Zeeman fields [14 were found to be 
governed by quantum laws due to quantum tunneling of 
magnetization [15]. In the experiments reported below, 
the spin reversal was triggered using a heat pulse at one 
end of the sample. 

Magnetic measurements were performed on three 
Mni2-Ac single crystals of dimensions ~ 0.3 x 0.3 x 2.1 
mm 3 , 0.35 x 0.35 x 1.75 mm 3 and 0.4 x 0.4 x 1.6 mm 3 
(samples A, B and C, respectively). The samples were 
mounted on a one-dimensional array of Hall sensors (ac- 
tive area 20 x 100 jim 2 with 200 /xm separation) (see Fig. 
[TJd). The three samples show essentially the same behav- 
ior; here we show data for sample C. Care was taken to 
align the sample and the Hall array (placed in the y-z 
plane) relative to each other and relative to the magnet 
axes, as shown in Fig. [I}}. The Hall sensors detect the 
stray field B x (see inset of Fig.j^J}), which is maximum 
at the reversal front. Measurements were taken at 0.4 
K in a 3 He refrigerator in a 3D vector superconducting 
magnet capable of producing bipolar bias fields up to 
1 T and bipolar transverse fields up to 5 T. A 6 V, 30 
ms pulse was applied to a thin film heater (R « 1.32 kfl 
at 0.4 K) placed at one end of the sample to trigger spin 
reversal. A dc current of 20 /iA was supplied to the ar- 
ray of hall sensors and the signal from each sensor was 
amplified by a factor of 1000, filtered, and recorded by a 
data acquisition card. 

At base temperature (0.4 K), the crystal was prepared 
in a fully magnetized state so that all the spins were 
aligned in one direction. A fixed H± was applied and H z 
was swept to the desired field and held. At this low tem- 
perature the spins are essentially blocked by the strong 
magnetic anisotropy of Mni2-Ac, and relax slowly to- 
ward equilibrium. In this non-equilibrium state, a volt- 
age pulse was then supplied to the heater to increase the 
temperature at one end of the crystal to initiate spin 
reversal. The same pulse was used to trigger the data 
acquisition card to record the magnetization signals of 
different Hall sensors during 1 s. The procedure was re- 
peated for several values of H z for a particular H±. 

Figure[2] shows the time evolution of the Hall sensor 
signals when a heat pulse is applied at t = to sample 
C with the magnetization prepared as described above. 
Two cases are considered: a zero bias field (H z = T, 
Fig.[2}i), and a large bias field (H z = 0.61 T, Fig.[2]o) 
for the same H± = 1.5 T. A propagating spin reversal 
region is observed in both cases; the change in magne- 
tization always starts at the heater end and propagates 
towards the other end (from 1 st to 6 th sensor, Fig. lip). 
As shown in Fig.[2j a peak in the signal corresponding 
to the enhancement of the local stray field (B x ) travels 
in the z-direction as the spin reversal propagates along 
the easy axis, ultimately reversing the magnetization of 
the entire crystal. The speed of propagation of the re- 
versing spins can be determined from the time difference 
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FIG. 2: Signals of the Hall sensor arrays as a function of time, 
a) At zero bias field, the peaks are very broad, the propagation 
time is long and the speed is not constant, b) At H z =0.61 
T, the peaks are sharp and spin reversal propagates at higher 
and constant speed. The inset shows stray fields from the 
sample when the reversal front is in the middle of the crystal. 
The maximum B x is at the position of the front. 

between pulses sensed by adjacent Hall sensors and their 
separation. 

While the curves in Fig.|2^i and Fig.[2]3 appear simi- 
lar, closer inspection shows that the reversal process at 
low and high bias are distinctly different. At zero (or 
low) bias, the time interval between adjacent peaks is 
relatively long and increases as the spin-reversal propa- 
gates, indicating that the spin reversal is slowing down; 
the time to completion is on the order of 80 ms; and the 
peaks are fairly broad and broaden further as the reversal 
proceeds. Here the spin reversal propagates slowly and 
its speed decreases because the magnetization relaxes to- 
ward equilibrium following the heat pulse; its progress 
through the crystal is governed by thermal diffusivity. 

By contrast, in the case of high bias shown in Fig.[2]3, 
the time between adjacent peaks is constant and on the 
order of 100 /is, corresponding to a constant speed of 
v ~ 2 m/s, consistent with the data measured by Suzuki 
et al. [5 ; the total time for reversal of the full mag- 
netization of the crystal is of the order of 1 ms; and 
the peaks in Fig.|2]3 are substantially sharper and do not 
broaden as they travel. The constant speed in the case of 
a large H z is a signature of magnetic deflagration, a self- 
sustained process driven by the Zeeman energy released 
by the spins as they reverse. 

We will now show that the crossover between these two 
regimes is surprisingly sharp. The speed of propagation 
of the spin reversal fronts was determined from the time 
interval between maxima of the Hall sensors in the middle 
of the crystal (sensor 3 and 4) and their spatial separa- 
tion. Note that while the same speed is obtained using 
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FIG. 3: a) The speed of propagation of the spin reversal as 
a function of H z for different H±. b) A blow-up of the data 
for small H z . The speed is small and nearly constant at very 
small H z and increases abruptly when a deflagration develops. 



Here T is the temperature, k is the thermal diffusivity, C 
is the heat capacity and Y = r exp [—U(h±,h z )/(kBT)] 
is the relaxation rate. For fields that are much smaller 
than the anisotropy field (hi = gjiBHi/2DS <C 1) [19] : 



U(h ± ,h z )^DS 2 (l-h z ) 2 



2h 



(1-^) 1/2 



. (3) 



For some applied fields and initial conditions, the ther- 
mal diffusivity cannot compensate for the increase in 
temperature due to the Zeeman energy released by the 
reversing spins; the sample temperature then rapidly in- 
creases and a magnetic deflagration develops. To study 
the ignition of magnetic deflagration, we consider a sys- 
tem of spins blocked in a given metastable magnetic 
state, increase its temperature (e.g., from 0.5 to 5K), 
and study the resulting nucleation process. Let us con- 
sider here for simplicity the temperature evolution in a 
nucleation volume, independent of coordinates, 



any pair of Hall sensors for large H Z1 the speed varies 
from point to point at low H Z1 and therefore between 
different Hall sensor pairs. As shown in Fig. [3] the speed 
of the spin reversal front changes abruptly at H z that 
depends on the transverse field. For a given H±, below 
a crossover bias field H co (for e.g., H co = 0.25 T for H± 
= T), the speed of propagation of the reversing spins 
is nearly independent of the bias field. However, above 
the crossover, the speed increases suddenly and depends 
strongly on the bias field. Figure^, a magnification of 
Fig. (3^i where the vertical scale has been expanded by a 
factor of 100, demonstrates this even more clearly. 

The peaks in the speed at H z ~ 0.45 T shown in 
Fig. |3^i are a clear manifestation of the important influ- 
ence of quantum mechanics on the dynamics of the sys- 
tem. These maxima are due to quantum tunneling at 
the resonant fields p~5j [16] whenever spin states on oppo- 
site sides of the anisotropy barrier have nearly the same 
energy. As noted earlier, tunneling enhances relaxation 
and effectively reduces the anisotropy barrier. As we in- 
crease H± we increase the tunnel splitting between levels, 
which results in a broadening of the resonance steps. A 
small peak can also be seen in Fig.[3]3 at the zero bias 
resonance. Note that the zero-field resonance is slightly 
shifted (~ 40 mT) due to internal dipolar fields [171 US] ■ 

To better understand the crossover between the two 
regimes described above, we develop a simple model 
based on previous theory [9j [10] . The diffusion-reaction 
system that describes the time evolution of the magne- 
tization, m, towards equilibrium, m eg , entails two pro- 
cesses: the Zeeman energy (i.e. heat) released into the 
sample by the spins as they reverse and the heat diffusion 
through the sample. The dynamical system of nonlinear 
partial differential equations reads: 



f = rh AE/C - 2k/R 2 (T - T ) 



(4) 



where T is the temperature of the volume under study, 
Tq is the temperature of the rest of the sample and 2R is 
the characteristic size of the nucleation volume. The neg- 
ative term in Eq. [4] is linear with temperature while the 
reaction term increases exponentially with temperature 
(fa oc exp [—[//(fee T)]); there is thus a competition be- 
tween these two terms and stationary solutions (T — 0) 
correspond to temperatures when the dissipation equals 
the released energy. 

As shown in the inset of Fig. [4] some magnetic field 
configurations (H = (H±,H Z )) corresponding to small 
bias fields have stationary solutions (red dotted curve). 
As we vary H (i.e., we vary the relaxation rate, m), 
the stationary points move closer to each other and 
eventually merge (green dashed curve) and disappear. 
When there is no stationary solution (blue solid line), 
the temperature derivative is always positive leading to 
an instability-a magnetic deflagration may develop. The 
condition for the crystal to lose stability at given T and 
H is T > and its crossover points are found by solving 



f = 0, df/dT = 0. 



(5) 



m 



-T(m 



f = rh AE/C + V • ftVT. (2) 



At a given H±, for H z < H co1 the energy released dur- 
ing the spin reversing process is small (or even zero in 
the case of zero bias). The dynamics of the propagating 
front observed in this case (Fig. [2^i) are determined by 
the diffusion of the heat supplied by the heater pulse; this 
corresponds to a thermal regime. Here the reaction term 
in Eq. [4] is negligible and the equation has a stable solu- 
tion at a temperature close to the sample temperature. 
Neglecting the reaction term in the full set of Eqs.[2j it is 
clear that the heat diffuses and the magnetization locally 
follows the sample temperature. 
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FIG. 4: Main panel: The boundary curves separating thermal 
and Zeeman regimes for various n. The blue dots are H co , 
derived from the data in Fig. [3] The inset plots T as a function 
of T for three different H. The solid and open black circles 
indicate stable and unstable solutions for the (1, 0.05) T field 
case. 

For H z > H co1 the energy released by the reversing 
spins is large, the reaction term dominates over the dif- 
fusion term and Eq. [4] displays an instability. In this Zee- 
man regime, the released energy drives the spin relax- 
ation to a deflagration with steady propagation, where 
the speed of the propagating front of reversing spins 
strongly depends on both H z and H±. 

The sharp transition in propagation speed from ther- 
mal to Zeeman regimes as a function of the applied mag- 
netic field is given by the condition in Eqs.[5] H± varies 
the activation barrier and does not change the released 
energy; it affects only m in the reaction term. H z affects 
both m and AE. 

In the main panel of Fig. [4] we plot the boundaries that 
separate the thermal regime from the Zeeman regime de- 
rived from Eq.[5] together with the experimental points 
(blue dots) obtained from Fig.[3j The different curves 
denote different initial metastable spin population (n = 
(M s — m)/(2M s ), where M s is saturation magnetization) 
that determines the initial state [26 . For fields above 
the boundary curve, the reaction term dominates and 
the crystal develops a magnetic deflagration; below the 
curves, thermal diffusion dominates and the spins evolve 
toward equilibrium through slow thermal relaxation. For 
the calculations we have bounded 2R by the smallest 
sample dimension of 0.4 mm, the heat capacity was taken 
to be C/ks = 1 (i.e. C = 8.3 J/(mol K)), close to a 
measured value [20], and the thermal diffusivity was es- 
timated to be k « 1 x 10 -4 m 2 /s [27], which agrees well 
with previously used values [21] . The free parameter is 
the temperature, To, to which the sample is heated by 
the voltage pulse; the best fit is obtained for To = 5 K. 

The smoothness and reproducibility of the data in 
Fig. [3] suggest that the appearance of instabilities leading 
to a magnetic deflagration when raising the temperature 



of the Mni 2 -Ac crystal is a well controlled process. Once 
the crystal is in the Zeeman regime the speed of the mag- 
netic deflagration can be controlled by the strength of 
the applied field. Additionally, equal propagation speeds 
are obtained for different combinations of applied fields, 
as seen in Fig. [3] where a particularly speed can be seen 
to correspond to different magnetic field configurations. 
Larger bias fields would produce larger amounts of heat 
released and higher temperatures for the overall process; 
while, large transverse fields that speed up the process 
by the same amount would heat the sample much less. 

We have measured the dynamics of spin reversal fronts 
in Mni2-Ac crystals in the presence of both bias and 
transverse magnetic fields. We have shown how these ap- 
plied fields can be used to control an instability that sepa- 
rates slow magnetic relaxation from rapid, self-sustained 
magnetic deflagration. 

We have also expanded the range of conditions under 
which magnetic deflagration has been observed, particu- 
larly to the case of small bias fields, which corresponds to 
small energy release in the deflagration process. This is a 
particularly interesting limit in which instabilities in the 
front are predicted to occur, such as pattern formation 
and oscillations of the front position as it propagates [22] . 
In the presence of large transverse fields, where quantum 
tunneling dominates over thermal effects, the instabilities 
in the magnetization dynamics may result in supersonic 
fronts [23] and lead to a deflagration-to-detonation tran- 
sition [24] [25] . The present experiments thus provide op- 
portunities to study-in a nondestructive manner-a large 
variety of instabilities in magnetic systems. 
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